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Quasinormal ringing of acoustic black holes in Laval nozzles is discussed. The equation for sounds 
in a transonic flow is written into a Schrodinger-type equation with a potential barrier, and the 
quasinormal frequencies are calculated semianalytically. From the results of numerical simulations, 
it is shown that the quasinormal modes are actually excited when the transonic flow is formed or 
slightly perturbed, as well as in the real black hole case. In an actual experiment, however, the 
purely-outgoing boundary condition will not be satisfied at late times due to the wave reflection at 
the end of the apparatus, and a late-time ringing will be expressed as a superposition of "boxed" 
quasinormal modes. It is shown that the late-time ringing damps more slowly than the ordinary 
quasinormal ringing, while its central frequency is not greatly different from that of the ordinary 
one. Using this fact, an efficient way for experimentally detecting the quasinormal ringing of an 
acoustic black hole is discussed. 

PACS numbers: 04.70.-s 



I. INTRODUCTION 

When the geometry around a black hole is slightly 
perturbed, a characteristic ringdown wave is emitted. 
This kind of phenomenon is known as quasinormal ring- 
ing, which is expressed as a superposition of quasinor- 
mal modes [l| . The central frequencies and the damping 
times of quasinormal modes are determined by the ge- 
ometry around the black hole, and thus the gravitational 
quasinormal ringing of a black hole is expected to play 
the important role of connecting gravitational-wave ob- 
servation to astronomy. 

Although quasinormal modes are themselves linear 
perturbations, they are in many cases excited after non- 
linear evolution of a black hole, such as the black hole for- 
mation after the merger of binary neutron stars. There- 
fore, in order to fully understand the emission process of 
quasinormal ringing in such a situation, one has to resort 
to numerical relativity, which needs extremely powerful 
computational resources. 

Interestingly, an alternative way to study the physics 
of black holes has been proposed. In a transonic fluid 
flow, sound waves can propagate from the subsonic re- 
gion to the supersonic region, but cannot in the oppo- 
site direction. Therefore, the sonic point of transonic 
flow can be considered as the "event horizon" for sound 
waves, and the supersonic region as the "black hole re- 
gion". For this reason, a transonic fluid flow is some- 
times called an acoustic black hole 0, Furthermore, 
it is shown that the wave equation for sound waves in an 
inhomogeneous flow is precisely equivalent to the wave 
equation for a massless scalar field in a curved space- 
time [1]. This implies that an acoustic black hole also 
has the quasinormal modes, and that one can study the 
physics of quasinormal ringing without resorting to nu- 
merical relativity, ultimately in laboratories. Although 
the quasinormal modes of acoustic black holes have been 
studied in some recent works [1, H, 0] , it was still an open 
issue whether they are actually excited and, if so, how 



they can be excited. 

In this paper, we show some results of our numeri- 
cal simulations to prove that the quasinormal ringing of 
sound waves is actually emitted from acoustic black holes 
in some situations, just as a real black hole emits the 
quasinormal ringing of gravitational waves. For future 
experiments in laboratories, we treat one of the most ac- 
cessible models of acoustic black holes: a transonic flow 
of air in a Laval nozzle [7,, i8,] . A Laval nozzle is a wind 
tunnel that is narrow in the middle, and is widely used 
to accelerate a fluid supersonically. With this nozzle, we 
can easily materialize a stable transonic flow in a labora- 
tory, just by exerting a sufficiently large pressure differ- 
ence between both ends of the nozzle. In addition, it is 
known that the sonic point is always formed at the throat, 
where the nozzle radius becomes the smallest. Thus, we 
do not need to search for the position of the sonic point 
in experiments. 

A quasinormal mode is characterized by a com- 
plex frequency (called quasinormal frequency) Un {n = 
0, 1, 2, • • ■ ), or equivalently, a pair of the central frequency 
Re(ci;„) and the quality factor Qn = Re(ci;„)/2|Im(aj„)|. 
The quality factor is a quantity that is propotional to the 
number of cycles of the oscillation within the damping 
time. In experiments, a damping oscillation like a quasi- 
normal ringing is inevitably buried in noise within a few 
damping times. Therefore, for detecting the quasinormal 
ringing efficiently, it is important to design a Laval noz- 
zle which yields a quasinormal mode with alarge quality 
factor. 

This paper is organized as follows: In Sec. II, we derive 
the Schrodinger-type wave equation for sounds, which de- 
scribes the "curvature scattering" of sound waves in an 
acoustic black hole, and compute the quasinormal fre- 
quencies of transonic airflow in Laval nozzles with known 
semianalytic methods. In Sec. Ill, we show the results of 
our fully-nonlinear numerical simulations to prove that 
the quasinormal modes of are actually excited in two dif- 
ferent types of situations. In Sec. IV, we discuss how the 
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quasinormal ringing changes when any outgoing waves 
are reflected at the upstream end of the apparatus. Sec. V 
is devoted to our conclusion. 



II. QUASINORMAL MODES OF ACOUSTIC 
BLACK HOLES IN LAVAL NOZZLES 

A. The Schrodinger-type wave equation 

In what foUows, we assume that the radius of the Laval 
nozzle r{x) varies sufficiently slowly along the axis of the 
nozzle (cc-direction) and the airflow in the nozzle can be 
considered as quasi-one-dimensional. In addition, we ne- 
glect viscosity and heat conduction throughout. 

The basic equations for quasi-one-dimensional flow of 
a perfect fluid are 



dtipA) + d^ipvA) = , 
dt{pvA)+d^[{pv^ +p)A]=0, 
dtieA)+d,[{e+p)vA]^0, 



(1) 
(2) 
(3) 



where p is the density, v is the velocity, p is the pressure, 
A is the cross section of the nozzle, and 



1 2 , P 
e= -pv + 

2 7 ^ f 



(4) 



is the energy density and 7 = 1.4 is the heat capac- 
ity ratio for air. Eqs. H]), and ^ are the equation 
of continuity, the momentum equation, and the energy 
equation, respectively. Insted of Eq. we can use Eu- 
ler's equation 



pidt + vdx)v = -d,;p, 



(5) 



which is derived from Eqs. (H]) and 

If the flow in the region under consideration has no 
entropy discontinuities (shock waves, contact surfaces, 
etc.), Eq. ^ can be replaced by the isentropic relation 



p (X p 



,1 



(6) 



Using this relation, Eq. ([5]) is reduced to Bernoulli's equa- 
tion 



5t$ + -(5,$)2 + /i(p) = 0, 



(7) 



where h{p) — J p^^dp is the specific enthalpy and $ 
J V dx is the velocity potential. 

Now we derive the linearized wave equation for sounds. 
To begin with, we replace (p, $) with the sum of the 
background part {p, $) and the perturbation part {Sp, 0): 

p = p + Sp, p>|<5/9|, (8) 
$ = $ + I^.l-I » |9,<^| , (9) 



Substituting Eqs. ^ and dH) into Eqs. ((T|) and 
and eliminating 6p, we obtain the wave equation for the 
velocity potential perturbation (p, 



- , (10) 



where a = pA^ v — dx^, and — \J dp /dp = \fyp{fi)Jp- 
Note that the sound speed depends on the local back- 
ground state, and is therefore a function of x. For 
A = const., Eq. (llOp is reduced to the one-dimensional 
version of Unruh's wave equation In the following, 
we omit the bars over the background quantities. 

Next, we consider the stationary background and 
Fourier transform (/)(i, a;) in terms of t: 



(j){t,x)e''^^dt. 



(11) 



Then, Eq. pH)) gives the time-independent diffcrntial 
equation for 4>^^, 



dx'' 



1 



(7(1 - Af2) 
1 



da; ( ^ 



ct(1 - Af2) 



r. d , 
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C2J 



1 - AI') ] + 2iuj^- 



= 0, 



d_ 

dx 



(12) 



where M = |u|/cs is the Mach number. In order to 
transform Eq. (|12p into the Schrodinger-type equation, 
we introduce a new function [x) and the "tortoise co- 
ordinate" a;* by 



H^{x)=g^'^ I dt e''^[*--^(")l0(t,a;), 
dx 



c.(l-M2)' 



where 



(J 

Cs 



\v\ dx 



(13) 
(14) 

(15) 
(16) 



and Cso is the stagnation sound speed, constant over the 
isentropic region. Then, Eq. (|12p can be rewritten as 



Vix*) 



0, 



where 



UJ 



K 



V{x* 
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9 d^g _ ^( dg^'^ 
2 dx*2 4 V dx* 



(17) 

(18) 
(19) 



The effective potential V{x*) has the dimension of 
(length) ~2 and characterizes the "curvature scattering" 
of sound waves on the acoustic black hole. In a Laval 
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nozzle, all of the profiles of the nondimensional quan- 
tities (p/po, vIvqi M, etc.) for the transonic flow are 
uniquely determined by a function A(x*)/Ath, where A,h 
is the cross-sectional area at the throat of the nozzle [3] . 
Hence, the structure of the effective potential V{x*) is 
also completely determined by A{x*)/At\i only. 

In this study, we consider a family of Laval nozzles 
which have the following form: 

A{x)=Tir{xf , (20) 
r{x) = roe - roo(l - /3) exp[-(x/2L)2"] , (21) 

where a is a positive integer and (3 = r(0)/roo £ (0,1). 
This type of Laval nozzle has a throat at a; = and 
two, infinitely long "tanks" with a constant radius r^c at 
X ^ L and x ^ —L . Note that the values of roc and 
L do not affect the structure of the nondimensionalized 
potential V{x*)L'^, and may be taken so that the quasi- 
one-dimension approximation [dr/dx <C 1) holds. Fig.[T] 
shows four examples of the form of the nozzle ±r{x) and 
the potential barrier V{x*) for the stationary transonic 
flow. 



and the curvature of the potential, respectively. Eq. (fT7|) 
with V{x*) = VpTix*) has a solution H^(x*) oc h^{x*) 
which behaves at x* — > ±oo as 



~ A(K)e-*''(^*-^o) + B(K)e^''(^*-^°) , x* 

where A{hi) and B{k) are defined by 

r(l - iK/K)T{-iHi/K) 



B{k) = 



r(s+ - iK/K)T{s^ - in/K) 
V{l-iK/K)T{iK/K) 



r(s+)r(s_) 



1 



- 1 



1/2. 



— oo , 
(25) 

-l-oo , 
(26) 



(27) 
(28) 

(29) 



Comparing Eqs. ^ and ^ with Eqs. ^ and ([26]). 
huj{x*) is found to be the quasinormal mode of the 
Poschl- Teller potential if A{k) = 0, or 



B. Evaluation of quasinormal frequencies 

Having obtained the form of the potential V{x*), we 
can compute the quasinormal frequencies by solving the 
Schrodinger-type wave equation (|17|) with the outgoing 
boundary conditions. 



X 

x 



-oo , 
-oo . 



(22) 
(23) 



In the case of a real black hole, the functional form of 
the effective potential is known analytically, and thus ap- 
plicable is Leaver's method which expands the wave 
function into a certain series and solves the continued 
fraction equation for the expansion coefficents. In the 
case of our acoustic black hole, however, the form of the 
potential can be known only numerically, and thus we 
have to apply an another method. 

In this study, we adopt the Poschl- Teller potential ap- 
proximation [ll| and the WKB approximation. 



1. Poschl-Teller potential approximation 

This method approximates a potential barrier with the 
Poschl-Teller potential, for which the scattering prob- 
lem can be exactly solved. The Poschl-Teller potential 
VpT(a;*) is defined by 



cos\i[K{x* - xl)] 



(24) 



where Xq is the position of the potential peak, and Vo = 
Vvt{xI) > and -2K^ = V^'t.{x^)/Vo < are the height 



(^0 



iK{s^ + n) 



iK 



(30) 



where n is a non-negative integer. With this relation, 
we can evaluate the quasinormal frequencies of a general 
potential V(x*) by approximating the form of the po- 
tential around with the Poschl-Teller potential, 
i.e., by setting = V{x*) and -2K^ = V"{xl)lV(xl) 
in Eq. pop . In general, this method gives an acurrate 
quasinormal frequency for small n, since such a mode 
is basically determined by the structure of V{x*^ at the 



maximum. 



2. WKB approximation 

This method constructs the quasinormal mode solu- 
tions by approximating H^{x*) with WKB functions in 
both sides of the potential barrier, matching them across 
the potential peak, and imposing the outgoing bound- 
ary conditions on it. The first-order WKB formula for 
quasinormal frequencies was derived by Schutz and Will 
[12| . and the higher-order versions have been developed 
by Iyer and Will flS*! , and Konoplya [ll| . 

The WKB formula for k„ reads 



Kl=V{x*)-^{n+-)[-2V"{x*)]i 

+ (higher-order correction terms) 



(31) 



where the higher-order correction terms depend on 
V'^^^\x^) and n (see Ref. [H 13). This method also 
gives an acurrate quasinormal frequency for small n. 
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FIG. 1: Forms of different Laval nozzles ±r{x) (thin curves) and the potentials V{x*) (thick curves) for stationary transonic 
flow, as functions of x. The fluid is assumed to flow in the negative x-direction, and thus 1 < is the supersonic region. Also 
indicated are the maximum point and the height of each V{x*) by the vertical and the horizontal dotted line, respectively. 
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TABLE I: Frequencies luq and quality factors Qo — 
Re(tt;o)/2Im(ti;o) of the least-damped quasinormal modes for 
different Laval nozzles, evaluated with the Poschl- Teller ap- 
proximation method (PT) and the WKB method up to the 
third-order (WKBl-3). The WKB results for a = 1 are not 
listed, since they do not converge as the WKB order increases. 



3. Results 

Table U shows the least-damped quasinormal frequen- 
cies llIq — CsqKq for Laval nozzles with different (a, /3), ob- 
tained by the Poschl- Teller potential approximation and 
the WKB approximation up to the third WKB order. 
For a > 2, the higher-order correction terms in Eq. (pij) 
converges well, and the Poschl- Teller values agree with 
the third order WKB values within ~ 10% accuracy. For 



a — 1, on the other hand, the correction terms are found 
to diverge as the WKB order increases, meaning that the 
WKB approximation breaks down. There is a possibility 
that the classical turning points across which the WKB 
functions are mathced are too distant from each other 
for these nozzles. To avoid confusion, we do not list the 
WKB values for a = 1 in the Table. 

We find that Re(Li;o)7 Ini(^o) ~ Cso/L, meaning that 
the wavelength and the damping length of the least 
damped mode are both determined by the curvature ra- 
dius L) of the Laval nozzle. We also find that the qual- 
ity factor Qo = Re(a;o)/|2Im(wo)| of the least-damped 
quasinormal mode is 0.5 <, Qq ^ 3 ^ 10, which implies 
that the quasinormal ringing damps out for only a few 
oscillation cycles. 

C. Numerical Simulations 

1. Settings 

In the numerical simulations, a Laval nozzle is assumed 
to have a form described by Eqs. and (HI]), and the 
flow in the nozzle is treated as quasi-one-dimensional. 
In order to solve Eqs. H]), ([2]), and ([3]) numerically, we 
adopted the MUSCL-Hancock scheme flSl. the second- 
order version of the Godunov scheme jl6l. [T7| . This 
scheme enables us to resolve shock waves generated in 
the flow in a good accuracy. The computational domain 
is set to —7L < X < 7L, divided into 800-1000 computing 
cells. 

In the case of an astrophysical black hole, quasinormal 
ringing appears when it is formed or when a test particle 
falls into it. As their analogues, we consider two types of 
situation for our acoustic black holes: (l) acoustic black 
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FIG. 2: Spacetime distribution of v (left) and log \v\ (right) 
obtained from a numerical simulation of type I. The nozzle 
parameters are set to (a,/3) = (2,0.01). The dashed curve 
and the dotted curve represent the position of the sonic point 
and the steady shock, respectively. 



hole formation and (11) weak- shock infall. In simulations 
of type I, the initial state is set to be stagnant and ho- 
mogeneous throughout the fluid, {p{x),v{x),p{x))t=o = 
(/90)0, Po)- At t — 0, the fluid in the nozzle begins to 
be "pumped out" from one end, at a rate high enough 
for the transonic flow to be formed. Fluid is allowed to 
enter freely from the other end so that the transonic flow 
can be stationary. In simulations of type II, on the other 
hand, flow is stationary and transonic from the first, and 
a weak shock wave is injected from the upper end to dis- 
turb the flow in the nozzle slightly. More infomation on 
the setting is summarized in Table HIl 
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FIG. 3: Numerical waveform Wobs(i) = v{t, 2.01/) (solid curve) 
obtained from the same simulation as that shown in Fig. [2] 
compared to the least-damped quasinormal mode (dotted 
curve) with the frequency ujo — (3.60 — 1.28j) Cso/L. 



the tanks, and thus the backreaction to the expansion of 
the fluid is much greater. The sonic point is formed at 
t ~ 13L/cso, and then approaches to the throat {x — 0) 
gradually. The discontinuity appearing in the supersonic 
region is the wavefront of a steady shock, which is pecu- 
liar to the stationary transonic flow in a Laval nozzle 0] . 
In the right panel of Fig. [21 we can see a ringdown wave 
emittied for t > 8L/csq from a; ~ L, around which the 
effective potential barrier is to be formed [see Fig.[TJd)]. 
This implies that the quasinormal modes are actually ex- 
cited during the transonic flow formation. 

In Fig. [31 the numerical waveform of v observed at 
X = 2.0L is plotted. The bumps at t ~ lOL/cso and 
13L/cso are due to the passing of shock waves generated 
in the throat region. For t <; ISL/c^o, we find that the 
waveform agrees very well with the least-damped (n = 0) 
quasinormal mode = 3.60 — 1.28z calculated by the 
Poschl- Teller potential approximation. From this fact, 
the ringdown wave appearing in Fig. [2] is proved to be 
the quasinormal ringing. 

The quasinormal ringing of transonic flow can be ob- 
served more distinctly in type II simulations, one of 
whose results is shown in Figs. [H and \5[ We find that 
the observed waveform is again in good agreement with 
the least-damped quasinormal mode. 



2. Results 

Fig.Hlshows the spacetime distribution of v and log \v\ 
obtained from the type I simulation for (a, /?) = (2, 0.01). 
At each position x, the fluid begins to move towards the 
lower boundary at t ~ (a: -|- 7L)/csq. In the "tank" re- 
gions (|a;| ^ L) the flow quickly relaxes to the stationary 
state, while in the "throat" region (|a;| ^ L) the flow set- 
tles down gradually, ejecting shock waves towards both 
tank regions. This difference reflects the difference in the 
cross section between the throat and the tanks; i.e., the 
cross section of the throat is much smaller than that of 



III. BOXED QUASINORMAL MODES: EFFECT 
OF THE WAVE REFLECTION AT A BOUNDARY 

Quasinormal modes are essentially characterized by 
the outgoing boundary conditions, Eqs. (l22l) and (I23p . 
In our acoustic black holes, Eq. (|^^ is always satisfied 
due to the presence of the acoustic horizon. On the other 
hand, whether Eq. (|23p is satisfied or not depends on 
the structure of the Laval nozzle in the upstream region. 
In the simulations shown in the last Section, for exam- 
ple, the upstream tank is assumed to be infinitely long, 
and therefore no wave is allowed to return to the poten- 
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type 


initial state 


boundary condition 




final state 


downstream [x — —7L) 


upstream {x = 7L) 


I 


{p,v,p) = (po, 0, po) 


(p,v,p) = (l.Opo, -O.lcso, 0.7po) 


transmissive 


stationary transonic 


11 


stationary transonic 


transmissivc 


transmissivc 


stationary transonic 



TABLE II: Settings of type I and type II simulations, 
of the flow, respectively. 



Po, Cso, and po denote the stagnation density, sound speed, and pressure 
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FIG. 4: Spacetime distribution of v (left) and log \v\ (right) 
obtained from a numerical simulation of type II, with nozzle 
parameters {a, (3) — (2,0.01). 
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tial barrier. In a such case, Eq. ([23|) is always satisfied 
and the potential barrier will always emit quasinormal 
ringing. In a real experimental configuration, however, 
the upstream tank has a finite length, and any outgoing 
wave emitted from the potential barrier is (perfectly or 
partially) reflected at the boundary, i.e., at the end of 
the tank. In such a case, Eq. ((23|) is satisfied only at 
early times, that is, before the wavefront of the reflected 
wave arrives at the potential barrier. Hence, in real ex- 
periments, the sound wave observed at late times will be 
different from an ordinary quasinormal ringing. 

Now we consider how the wave reflection at the bound- 
ary alters the quasinormal ringing in the Laval nozzle. In 
the following, we assume that the upstream tank has a 
reflective boundary wall perpendicular to the axis of the 
nozzle at X = Xc ^ xq. Then, the appropriate boundary 
conditions for are 



~oo, (32) 
<, (33) 



where TZ^j G [0, 1] is the reflection coefficient of the 
boundary wall. The boundary wall is perfectly reflective 
{v\x* cx; dxHi^lx* = 0) for 7?,^ = 1, and partially absorb- 
ing for < 7^ < 1. The solutions of the wave equa- 
tion ()17p with these boundary conditions are no more 
quasinormal modes in a narrow sense, but are yet ex- 
pected to be "quasinormal" , with complex characteris- 
tic frequencies. For this reason, we shall refer them to 
"boxed" quasinormal modes after Refs. [sL [l8l|. 



A. Evaluation of boxed quasinormal frequencies 

In order to calculate the frequencies of boxed quasinor- 
mal modes, we approximate V{x*) with the Poschl- Teller 
potential used in Sec II. Comparing the upstream bound- 
ary condition (|33p with the asymptotic form (|26p of the 
particular solution h^^^x*), we obtain the characteristic 
equation for k, 



r{s+)r{s,)r{-iK/K) 



r(s+ - iK/K)r{s- - iK/K)T{iK/K) 



2iKA 



(34) 



where A = x* — a:g is the distance between the boundary 
wall and the peak of the potential. For simplicity, we 
assume the reflection coefficent TZ^ is independent of u; 
and rewrite it as TZ^^ = TZ. We find a number of solutions 
for fixed A, and the number is considered to be infinite 
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FIG. 6: Boxed quasinormal frequencies for (a, l3) = (2, 0.01) 
and TZ = 0.089, as functions of Xc- The curves indicate eight 
of the numerical solutions to Eq. (I34p . ujh,m{m — 0, • • ■ ,7). 
For comparison, the least-damped ordinary quasinormal fre- 
quency LOo is indicated by the thick horizontal line. The sym- 
bols 'o', 'x', 'A', '<>', and represent the "best fitting" fre- 
quencies ojinin for the waveforms obtained through type F 
simulations (see Sec. III.B). 
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FIG. 7: Numerical waveform Uobs(t) = v{t,2.0L) (solid) for 
an early time (upper panel) and a late time (lower panel), 
obtained from the type F simulation with (q, /3, Xc, 7?.) = 
(2, 0.01, 4. OL, 0.089). For comparison, the least-damped qua- 
sionormal mode (dotted), a superposition of three boxed 
quasinormal modes (dashed), and the "best fitting" function 
(dot-dashed) are also plotted. 



for a reason mentioned below. We refer to the solutions 
as Kb,m, where m = 0, 1, 2, • ■ • are labeled in the order of 
increasing Re(Kb,m), not of increasing |Ini(Kb,m)|- 

In Fig. [SI we plot some of the numerical solutions of 
Eq. dMl), ujm (m = 0, ■•■ ,7), for = (2,0.01) and 

n = 0.089 as functions of Xc A + 0.63L). We find 
that the central frequencies Re(ci;b,m) approximately sat- 
isfy the relation Re(a;b,m+i) — R.e(wb,m) — t^Cso/ A, which 
is the same as the one satisfied by the normal modes 
in a cavity with length A. This means that the boxed 
quasinormal modes are basically the same as the normal 
modes in a cavity; the only difference is that the "cavity" 
is leaky, i.e., the boundary wall and the potential barrier 
transmit sound waves partially. For this reason, the num- 
ber of the boxed quasinormal modes is considered to be 
infinite as is for the normal modes. 

We also find that the boxed quasinormal modes with 
Re(wb,m) ^ R-e(wo) have the damping rates |Im(ti;b,m)| 
that are much smaller than that of the least-damped 
quasinormal modes, |Im(a;o)|- From this fact, the ringing 
observed at late times is expected to damp much more 
slowly than the quasinormal ringing at early times. In 
fact, it is found from the results of the numerical simu- 
lations that the modes with Re(u;b,m) — Re(u;o) are the 
most strongly excited, and as a result the late-time ring- 



ing has the quality factor much smaller than that of the 
early-time quasinormal ringing (see Sec. III.B). We also 
calculated the frequencies for to > 8 and checked that 
their damping rates are higher than that for m — 7 (at 
least in the range 2L < Xc < 6L). Therefore, the modes 
with large to will not dominate the late-time ringing, even 
if the number of such modes is infinite. 



B. Numerical simulations 

In order to verify that the late-time waveform is acually 
characterized by the boxed quasinormal modes, we per- 
form modified versions of type I simulations, to which we 
refer as type F simulations. In type F simulations, the 
upstream boundary is set to be reflective (0 < 7?. < 1) 
and located at x = Xc(> xq), while the initial state in 
the nozzle and the condition at the downstream bound- 
ary are the same as those of type I. A summary on the 
setting of type I' simulations is shown in Table Hill 

Fig. [7] shows the numerical waveform of v at 
Xobs — 2.0L obtained from the type I' simulation with 
{a,i3,Xc,TZ) = (2, 0.01, 4.0L, 0.089). At early times {t 
lAL/cso), the waveform looks identical to that shown 
in Fig. [21 and is well fitted to the least-damped quasi- 
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type 


initial state 


boundary condition 
downstream [x — —7L) upstream [x = Xc) 


final state 


I 


ip,v,p) = (po, 0, po) 


{p,v,p) = (l.Opo, -O.lc^o, 0.7po) reflective (0 < 7^ < 1) 


stationary transonic 



TABLE III: Setting of type F simulations. TZ denotes the reflection coefficient at the upper boundary. 



normal mode again. In fact, t = lAL/csO is just the 
time at which the bump observed at i ~ lOL/cso re- 
turns to a; = Xobs = 2L after reflected at the bound- 
ary (x = Xc = 4L). At late times (t ^ ML/cso), the 
waveform is no longer fitted to the ordinary quasinormal 
modes. In order to verify that the late-time waveform is 
a superposition of the boxed quasinormal modes, we de- 
compose the waveform into damped sinusoids as follows: 
first, we calculate the function 



Eiuj,C,Vo 



*2 r -| 2 

Vohsit) - u{uj,C,Voo;t) dt, (35) 



tl 



u{u;,C,Vca;t) = Woo +Re 



^ (-.(i)g-ic^O)(t-tl) 



, (36) 



where a; = (w^^), • • • c = {C^^\- ■ ■ , C^^'— )), 

and Voo are fitting parameters, and ti, t2 are taken to be 
large enough. Then, we search for the set of parameters 
(w, C, Woo)min that minimizes E. In this example, setting 
Jmax = 5, tl — 19L/cso, and t2 — 35L/cso, we obtain 
Winin = (2.61 - 0.348i, 3.54 - 0.364i, 4.13 - O.U9i)cso/L. 
In the following, we refer to u(tt',nin, Cmin, Woo.min; ^) as 
the "best fitting" function. In fact, the waveform is 
well approximated by the "best fitting" function for 
t ^ 18L/Cso, as shown in the lower panel of Fig. [T] 
In addition, we find that u^min agree with three of the 
boxed quasinormal frequencies calculated from Ea. ([M)) . 
(wb,3,Wb,4,Wb,5) = (2.47 - 0.356i, 3.34 - 0.353i, 4.15 - 
0A30i)cso/L, within 6% accuracy. For comparison, a su- 
perposition of the three boxed quasinormal modes is also 
plotted in Fig. [71 It is noted that the quality factor of 
'^min is Q = 4.44 on average, while that of the least- 
damped ordinary quasinormal mode is Qq = 1.41 (in the 
Poschl- Teller approximation). It is quite surprising that 
a weakly-reflecting wall with the reflectance |7^p ^ 1% 
enhances the quality factor of the ringing by ~ 300%. 

We carry out type I' simulations for different Xc and 
obtain Wmin with the same analysis as above, whose re- 
lults are plotted in Fig. [6l In the analysis, we need to 
increase the number jmax of the damped sinusoids with 
increasing Xc, since the number of the boxed quasinor- 
mal modes with small |Im(a;b)l s-lso increases. It is found 
that the observed frequencies corresponds to those of the 
boxed quasinormal modes with Re(a;b,m) close to Re(a;o)j 
implying that the boxed quasinormal modes which are 
"close" (in the frequency space) to the least-damped or- 
dinary quasinormal mode are strongly excited. This is 
narural since the boxed quasinormal modes are excited 
by the quasinormal ringing at the early time. Roughly 
speaking, the late-time ringing inherits the central fre- 
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FIG. 8: Numerical waveform Vohsit) = v{t,2.0L) (solid) for 
an early time (upper panel) and a late time (lower panel), 
obtained from the type F simulation with (q, /3, TZ, Xc) = 
(2,0.01, 1.0,4.0L). 



quency from the quasinormal ringing at the early time, 
while its damping rate is greatly supressed due to the 
reflection at the boundary. This is the main reason why 
the quality factor (oc the central frequency divided by 
the damping rate) is greatly enhanced by the reflective 
boundary. 

Up to here, we have dealt with the effect of a weakly- 
reflecting boundary {TZ <^ 1) and shown that the our 
formula (|34p gives the boxed quasinormal frequencies in 
a good accuracy for that case. Now we show that the 
formula (IM)) is also valid for a perfectly-reflecting bound- 
ary (TZ =1). Fig. [7] shows the numerical waveform of v 
at Xobs = 2.0-L obtained from a type I' simulation with 
{a,(3,TZ,Xc) = (2,0.01, 1.0, 4.0L). We carry out for this 
waveform the same analysis as above, and obtain iVrnin — 
(2.59-0.00546i, 3.45-0.0430i, A.28-0.102i)cso/L. These 
values agree with three of the boxed quasinormal fre- 
quencies calculated from Eq. ([M)) . (wb,3, Wb,4, Wb.s) = 
(2.48 - 0.00457i, 3.36 - 0.0302i, 4.20 - O.iui)cso/L. The 
averaged quality factor of Wmin is Q = 99.4, which is far 
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greater than Qo- 

IV. CONCLUSION 

We have studied the quasinormal modes of acoustic 
black holes in Laval nozzles both semianalytically and 
numerically. We have found that the quasinormal modes 
are actually excited in response to perturbations to the 
transonic flow, but also found that the quality factor 
is much smaller than 10, i.e., the quasinormal ringing 
damps within only a few oscillation cycles. In an ac- 
tual experiment, however, the purely-outgoing boundary 
condition will not be satisfied due to the reflection at 
the end of the apparatus, and the observed ringing will 
be expressed as a superposition of "boxed" quasinormal 
modes, not of the ordinary quasinormal modes. The 
quality factor of the late-time ringing can exceed that of 
the least-damped ordinary quasinormal mode, since the 
central frequency is basically inherited from the early- 
time, ordinary quasinormal ringing while the damping 
rate is greatly suppressed by the reflective boundary wall. 
Therefore, if the inner wall at the uppstream end of the 
experimental apparatus is well arranged so that it reflects 
sound waves as perfectly as possible, the detection of the 
quasinormal ringing in a laboratory will be done quite 
efficiently. 

In this study, we have dealt with quasi-one-dimensional 



transonic flows only. More generally, it is possible to 
consider a flow with a nonvanishing azimuthal veloc- 
ity component, an acoustic analogue of a rotating black 
hole. It is known that the quality factor of the least- 
damped quasinormal mode is larger for rotating black 
holes than for non-rotating. Hence, the quality factor of 
the quasinormal modes for a "rotating" acoustic black 
hole (transonic flow with azimuthal rotation) may be 
large than that for "non-rotating" ones. Also, it is consid- 
ered that a reflective shell that encloses a rotating black 
hole can explode due to the super-radiant scattering by 
the black hole; such a system is well-known as a "black 
hole bomb" [2Q(| . Similarly, a rotating acoustic black hole 
with a reflective boundary wall is expected to behave as 
an "acoustic black hole bomb" [J. We are planning to 
perform numerical simulations for rotating acoustic black 
holes in Laval nozzles and investigate these issues in the 
future works. 
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